Probing quasiparticle excitations in a hybrid single electron transistor 
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We investigate the behavior of quasiparticles in a hybrid electron turnstile with the aim of improving its performance 
as a metrological current source. The device is used to directly probe the density of quasiparticles and monitor their 
relaxation into normal metal traps. We compare different trap geometries and reach quasiparticle densities below 3 jim"^ 
for pumping frequencies of 20 MHz. Our data show that quasiparticles are excited both by the device operation itself 
and by the electromagnetic environment of the sample. Our observations can be modelled on a quantitative level with 
a sequential tunneling model and a simple diffusion equation. 



Applications of superconductors generally rely on the fact 
that electronic excitations can be generated only if energy 
higher or equal to the gap energy A is available. Hence 
the number of the excitations is ideally exponentially small 
at low temperatures and the properties specific for super- 
conductors appear. If the number of excitations, typically 
characterized by the density of quasiparticles, increases, 
the superconducting features degrade. Such an effect has 
been studied in several devices such as superconducting 
qubits^"^, superconductor-insulator-normal metal-insulator- 
superconductor (SINIS) microcoolers^"^ and kinetic induc- 
tance detectors^' In this letter we focus on the effects of 
quasiparticles on a hybrid single-electron transistor. We use 
the Coulomb blockaded transistor for direct and simple prob- 
ing of the quasiparticle excitation density. The device operates 
as a charge pump and is a promising candidate for the realiza- 
tion of a metrological current source^ ^ We show experimen- 
tally that the quasiparticle excitations limit the current quan- 
tization. By optimizing the quasiparticle relaxation, however, 
we estimate that it is possible to reach metrological accuracy. 

In order to observe how quasiparticles influence the perfor- 
mance of the turnstile we designed samples (type A) where 
the quasiparticle relaxation in normal metal traps was pur- 
posefully delayed by extending the bare superconducting lines 
that connect the junctions to the traps. The beginning of this 
isolated superconducting line can be seen in the scanning elec- 
tron microscope image of the turnstile in Fig. 1. (a) where it 
connects to the normal metal island via oxide junctions that 
appear as lighter areas. Its extension is visible in (c) all the 
way through to the wide traps of overlapping superconductor 
and normal metal separated by an oxide layer. The sample 
shown in (b) (type B) has wide leads with normal metal traps 
close to the junctions to enable efficient quasiparticle evacua- 
tion. The samples were fabricated with the standard electron- 
beam lithography and shadow mask technique We com- 
pare the behaviour of quasiparticles in SINIS turnstile samples 
with different geometries. The length of the isolated supercon- 
ducting line (given by the separation of the transistor junction 
from the trap) was varied between / = 200 nm and / = 20 jim. 
Measurements were performed in a dilution refrigerator at a 
base temperature of approximately 60 mK. 

Figure 1 (d) shows the current- voltage characteristics of the 




FIG. 1 . (a) Scanning electron micrograph of the / = 20 jim type A 
SINIS20 turnstile showing the superconducting Al leads connected 
to the Cu metal island via AIO2 barrier junctions and the gate elec- 
trode used to regulate the potential of the island, (b) shows the type B 
SINISopen geometry sample where the normal metal traps are within 
200 nm of the junctions and the leads open up straight from the junc- 
tion. In (c) the full length / of the superconducting line of sample 
SINIS20 between the junction and the trap is visible as well as the 
quasiparticle traps formed by overlapping Cu and Al shadows. A 
sketch of the basic measurement circuitry is depicted and (d) shows 
the IV characteristics of the turnstile. 



turnstile sample presented in (a) and (c). The black line shows 
the current through the turnstile when the bias voltage is 
swept across the superconducting gap and the gate potential 
is varied between the gate open {Ug = 0.5) and the gate 
closed (rig = 0) states. Simulations based on sequential tun- 
nelling are fitted to the data in order to extract the parameters 
of the superconducting gap A = 216 jieV, the charging en- 
ergy of the island Ec = 0.74 A and the tunnelling resistance 
of the junction Rt = 91 kQ, specific to each sample (see Ta- 
ble 1 for parameters of all samples measured). These sim- 
ulations are shown in blue for the gate open and red for the 
gate closed state. On this coarse level, the heating of the su- 
perconducting leads does not make a significant contribution 
and can be neglected. The overheating of the normal metal- 
lic island (thickness 30 nm) is taken into account by consid- 
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FIG. 2. (a) and (b) show the normaUzed current through the oxide 
trap turnstile of geometry / = 20 jim smoothed over 5 data points as a 
function of Vb for gate voltages varying between = and rig = 0.5. 
The charge stability measurement plotted in (a) was performed in a 
sample stage with a single cover; (b) shows measurements of the 
same sample in an indium sealed two-cover stage. In (c) and (d) we 
plot the normalized current computed with a simulation based on se- 
quential tunneling for the situations of (a) and (b) respectively. The 
inset of (d) shows the quasiparticle density Ugp inferred from the tem- 
perature reproducing the measured current inside the superconductor 
gap as a function of the length / of the isolated superconducting line. 



ering the electron-phonon coupling with material parameter 
value of S = 2 • 10^ WK"^m"^ which is consistent with val- 
ues obtained in previous experiments In the subgap and 
pumping experiments we have the opposite situation: the is- 
land does not heat since the signals are much smaller but as 
we look at small differences, the overheating of the supercon- 
ductor starts to have an effect. 

We now investigate the current of the turnstile more closely 
for bias voltages within the superconducting gap. These mea- 
surements allow the direct probing of the quasiparticle den- 
sity. This is not possible with voltage biased NIS junctions 
that probe the excitations more indirectly^"^'^^ and correspond 
to the gate open case in these experiments. Figure 2 (a) shows 
the current as a function of Vb and Vg. With no quasipar- 
ticles present we would expect zero current for low biases. 
However, we observe a current pattern periodic in Vg within 
the gap with currents rising up to 10 fA. We find that we 
can reproduce the measured current pattern with a simulation 
shown in panel (c) using a high superconductor temperature 
Ts that gives rise to a quasiparticle population above the gap. 
We compute the average current through the left junction via 
I = e Y^ni^Liin) - TiL(n))P(n, t) where Tu{il) is sequential 
tunnelling rate to (from) the island through the left junction 
and Pin, t) is the probability of the system being in the charge 
state n of the island. The temperatures used to fit these data are 
Tg = 205 mK and the normal metal temperature = 92 mK. 
The current- voltage characteristics are surprising: At degener- 
acy (half integer rig) no net current flows, whereas in Coulomb 
blockade (integer rig) we obtain a finite current. The simula- 



tions give us an insight into the on-going processes: at degen- 
eracy the hot quasiparticle excitations lying at high energies 
are able to tunnel in both directions equally, hence there is no 
net current. In Coulomb blockade, the tunneling of a quasi- 
particle excitation is followed by a fast relaxation to the lowest 
lying charge state. The relaxation always happens in the for- 
ward direction given by and leads to a net current through 
the device. These features are a strong indication of having 
quasiparticle excitations as the source of the sub-gap current 
in the device. 

We measured the same sample in two different sample 
stages. Instead of having an enclosed stage with one metal- 
lic cover as in the measurement displayed in Fig. 2 (a) we 
also used an indium sealed double hermetic metallic sample 
stage. In the latter case, the radio frequancy (RF) line had 
an additional 20 cm thermocoax cable to enhance the sample 
shielding. This was not included in the wiring of the single 
cover stage. However, as the RF line is not electrically con- 
nected to the turnstile, we do not expect this to influence the 
direct heat conduction to the sample. The dc wiring of the two 
setups was similar, made of approximately two meters of ther- 
mocoax cable. The two sample stages are thermalized to the 
cryostat base temperature in an identical way and the sole pur- 
pose of this sealing is to create an RF shield to the sample. The 
result is shown in panel (b). In this case no sub-gap current 
can be resolved. This behavior was fitted with temperatures of 
< 167 mK and = 72 mK which we present in Fig. 2 (d). 
The inset in Fig. 2 (d) displays the quasiparticle density riqp 
inferred from the relation riqp = 

^[2n D(Ef)A ^fkBTj^ e'^^^^^ which is valid at low temper- 
atures kBTs « A for three different samples with varying 
distance to the trap and measured with the poorly filtered sam- 
ple stage. We use D{Ep) = 1.45 x lO'^^m"^/"^ as the normal 
state density of states at the Fermi energy . We observe a 
monotonous increase in riqp with increasing distance of the 
trap from the junction. We deduce that the presence of envi- 
ronmentally excited quasiparticles is determined by both the 
trap relaxation rate and the diffusion rate through the super- 
conducting line. 

Next we turn to the dynamic case of the turnstile opera- 
tion. In addition to environmental excitation, quasiparticles 
are now injected to the superconducting leads once in every 
pump cycle. The pumping frequency thus allows us to control 
the injected power and the number of quasiparticles. Figure 3 
shows current plateaus measured on type A SINIS20 sample. 
Three bias voltages were chosen around the optimum opera- 
tion voltage of eVh = A^^ ranging from 0.8A/^ to 1.6A/^. Two 
main effects were observed: a slight overshoot for the highest 
bias voltage at the beginning of the plateau and a spreading of 
the plateau value for different V^. Simulations including sin- 
gle electron and two-electron Andreev processes were used to 
model the process. The peak at the beginning of the plateau 
can be reproduced by allowing for an Andreev current with a 
conduction channel area of 30 nm^. This value is taken from 
previous experiments with similar samples The spread- 
ing of the plateaus in Vt is caused by quasiparticle excitations 
and follows the expected behaviour with respect to variations 
in geometry and pumping operation as described below. 
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FIG. 3. Current plateaus under pumping. Current through the turn- 
stile as a function of the gate modulation amplitude measured on 
the type A SINIS20 sample at (a) 1 MHz, (b) 10 MHz and (c) 20 
MHz sinusoidal modulation frequency and on the type B sample 
SINISopen in (d) at 20 MHz. The pump was operated at bias volt- 
age values of eVb = 0.8A (blue dots), eVt> = 1.2A (green dots) and 
eVb = 1.6A (red dots). Fits from simulations including two-electron 
Andreev processes are displayed as solid black lines, dashed black 
lines are simulations only including single electron processes. 



The spread increases with the pumping frequency / from 
A/ = 8 fA for 1 MHz (panel (a) of Fig. 3) to A/ = 15 fA for 
10 MHz (panel (b)). The injected power Pinj increases with 
/ and we model this increase in mean number of quasiparti- 
cles in the superconductor by raising its temperature Ts . The 
bias-dependence of current on the plateau is not influenced by 
Andreev currents and we can thus fit it simply by ascribing 
it to the quasiparticle number. The black lines correspond to 
simulations with quasiparticle densities of riqp =41.1 jim"^ at 
1 MHz and Uqp = 82.6 jim"^ at 10 MHz. 

To demonstrate that the quantization can be improved 
by enhancing the quasiparticle relaxation we measured the 
pumping plateaus of a sample with broader leads, see 
Fig. 1 (b). Panels (c) and (d) in Fig. 3 show the comparison 
of poor and good quasiparticle trapping respectively. The cur- 
rent quantization improves by two orders of magnitude when 
using leads with enhanced quasiparticle relaxation. 

To model the relaxation of the quasiparticle excitations we 
consider their diff'usion in a thin superconducting line. The 
heat diff'usion equation is 

V.(-^5Vr) = -Arap, (1) 

where the thermal conductivity of the superconductor is 
= ^ (^) e~^'^^^ i^{)T I pn with the Lorentz number Lq 
and the normal state resistivity p„^^'^^. The heat is re- 
moved from the superconductor by a normal metallic trap 
to which it is tunnel coupled. The trap removes the heat 
Arap = ^ Ens (E)(fN(E)-fs (E)) dE per unit area where 
(Tt is the conductance of the trap per unit area and d = 22 nm 
is the thickness of the superconducting film. The power Pinj 



injected into the line during turnstile operation sets a bound- 
ary condition to the beginning of the line: Pinj = A • (-ksVT), 
where A = wd is the cross-sectional area of the lead with 
width w. We rewrite Eq. (1) in terms of quasiparticle density 
riqp by considering only the strong exponential dependencies 
on T to obtain 

V^riqp = A'^iriqp - Hqpo), (2) 

where = ^^^^ (¥) ^ ^ip^ quasiparticle den- 
sity of the superconductor when it is fully thermalised to the 
normal metal. Next we solve Eq. (2) for three diff'erent sec- 
tions of the bias line for type A samples. In the first part we 
treat the bare aluminium line of length / and constant width 
wi with no quasiparticle trapping (ptrap = 0)- The next part 
deals with the widening of the lead from wi to W2 where we 
also neglect quasiparticle trapping since its contribution in this 
section of the lead is small for our samples. Then, in the last 
part, the lead continues with a constant cross section of W2 and 
is in contact with the quasiparticle trap: ptrap 9^ 0. As a result, 
we obtain the quasiparticle density at the junction to be 

Hqp = ;=== - / + wi log — +.!— . (3) 

The first term is the diff'usion in the bare aluminium wire, the 
second term is the spreading to the wider line and the last 
term arises from the relaxation to the trap. Similarly, we can 
solve the diff'usion for the type B sample with opening bias 
lines. Here we assume that the line starts at radius tq and 
the injected power is distributed evenly to all directions. We 
take ro = 70 nm so that the area the power is injected into, 
|ro<i = (50 nm)^, matches the junction area of the sample. 
The quasipartice density at the junction is then 

^ ^e^D{EF)pn Pinj Ko(ro/A) 
^/2Ak^ OrodKiiro/AY 

where Kn is the modified Bessel function of second kind and 
the opening angle of the line. 

We now compare the quasiparticle relaxation in diff'erent 
sample geometries. From fits to measurements similar to 
those shown in Fig. 3 we extract riqp as a function of /. These 
values are displayed as dots in Fig. 4 for the SINIS20 (blue), 
the SINIS5 (green) and the SINISopen samples (red). Us- 
ing the diff'usion model described above we calculate riqp as 
a function of / corresponding to an injection power on the 
plateau given by P^j = efA (solid lines). The injected heat 
calculated from the simulations deviated less than 10 % from 
this power even at the highest frequencies measured. We used 
p„ = 31 nQm as the normal state resistivity for all samples. In- 
dividual sample parameters used in the simulations are listed 
in Table 1 . The measured densities clearly show a linear be- 
havior with increasing / as predicted by the model (see Eqs. 
(3) and (4)). As indicated by the slopes for increasing length / 
of the isolated line, the further the quasiparticle trap lies from 
the oxide junction and the thinner the connecting line is, the 
slower the relaxation process is and the more the supercon- 
ductor is heated. For the samples SINIS20 and SINIS5 a fi- 
nite quasiparticle density is observed even in the absence of 
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FIG. 4. Quasiparticle density as a function of the driving gate voltage 
/. Fits to measurements including first and second order tunnelling 
processes are displayed as dots and quasiparticle densities derived 
from the diff'usion model are displayed as solid lines. We compare 
two type A samples with delayed relaxation, SINIS20 (blue) and SI- 
NIS5 (green) and one type B sample SINISopen (red). Diagrams 
of the sample geometries are displayed in the upper part of the fig- 
ure, dark colours representing quasiparticle traps, lighter areas the 
isolated superconducting lines and the small black bulks the normal 
metal islands. Below, a diagram shows the two pathways to quasi- 
particle excitation, during pumping via the normal metal island and 
through interactions with the electromagnetic environment. 



injected power. This density is of the same order of mag- 
nitude as the leakage currents shown in Fig. 2 and it points 
towards interactions with electromagnetic radiation from the 
environment. The SINISopen sample clearly shows the most 
efficient quasiparticle relaxation with only a few quasiparti- 
cles per [im"^ at 50 MHz driving frequency. 



TABLE I. Sample parameters 
Sample Edl^ Rt [kQi] / [jim] wi [nm] W2[!im] crj[(Qm^)~^] 



SINIS20 


0.74 


91 


21 


120 
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2.3x10^ 


SINIS5 


0.7 


81 


5.4 


120 


5 


2.7x10^ 


SINISl 


0.7 


85 


1.1 


160 


5 


9.5x10^ 


SINISopen 


1.2 


170 








l.OxlO'^ 



We have investigated the process of quasiparticle relaxation 
in hybrid SINIS turnstiles for various geometries. Both direct 
current and gate-driven pumping measurements can be well 
understood and simulated using models based on first and sec- 
ond order tunnelling processes. We were also able to model 
the relaxation of the excitations by a simple diflfusion model. 
We find that the main sources for quasiparticles in the single 
electron transistor are those injected above the superconduct- 
ing band gap during turnstile operation and those excited by 
radiation from a hot environment. In the best structure stud- 



ied in this Letter we reached an accuracy 61 1 1 of the order of 
1 • 10""^. We estimate that accuracy better than 10"^ at 50 MHz 
would be obtained with the following improvements: The alu- 
minum should be made an order of magnitude thicker and the 
trap ten times more transparent. Alternatively, the bias leads 
can be extended to the third dimension. The resistance of the 
junctions should be increased by a factor of four and charg- 
ing energy to Ec > 2A. Finally, the density of environmen- 
tally activated quasiparticles needs to be reduced to a level 
^qp ^0.1 yum"^ which has been demonstrated experimentally 
in a recent work^^. With these realistic modifications we ex- 
pect to reach metrological accuracy in turnstile operation. 
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